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We report on a search for extremely-high energy neutrinos with energies greater than 10 6 GeV 
using the data taken with the IceCube detector at the South Pole. The data was collected between 
April 2008 and May 2009 with the half completed IceCube array. The absence of signal candidate 
events in the sample of 333.5 days of livetime significantly improves model independent limit from 
previous searches and allows to place a limit on the diffuse flux of cosmic neutrinos with an E~ 2 spec- 
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I. INTRODUCTION 

Cosmogenic neutrinos, the daughter particles of the Greisen-Zatsepin-Kuzmin (GZK) process in which the highest 
energy cosmic-rays interacting with the cosmic- microwave background [H,[3|, may give a unique picture of the Universe 
in the highest energy regime. Cosmogenic neutrinos carry information about the sources of the highest energy cosmic- 
rays, such as their location, cosmological evolution, and cosmic-ray spectra at the sources. Various cosmogenic neutrino 
models which assume primary cosmic-ray protons predict neutrino fluxes Efy > 10~ 4 GeVm _2 sec _1 sr _1 in the 
energy range 10 8 GeV <E V < 10 10 GeV, which implies that the 4-7T solid angle averaged neutrino effective area divided 
by energy A v /E v must be larger than 1CT 5 m 2 /GeV (e.g. A v > 10 3 m 2 at 10 s GeV and A v > 10 4 m 2 at 10 9 GeV) to 
detect several cosmogenic neutrinos every year. 

Several techniques have been used to realize such huge detection volumes for these extremely-high energy (EHE) 
neutrinos. Air-shower detectors search for neutrino induced young inclined showers Q or Earth-skimming events 
initiated by tau neutrinos Q. Radio Cherenkov neutrino detectors search for radio Askar'yan pulses in a dielec- 
tric medium as the EHE neutrino signature [9|-llll]. Underground neutrino telescopes, such as IceCube, deployed in 
transparent naturally occurring media [HI, [TJ can detect EHE neutrino interactions through the strong Cherenkov 
radiation emitted by the charged secondary particles. This technique is well established for observations of astro- 
physical neutrinos in the MeV to GeV energy region [3 EH , and can also be utilized to search for cosmogenic EHE 
neutrinos with an appropriate background rejection method. In a neutrino telescope, an EHE neutrino interaction 
would be identified by the extremely high number of Cherenkov photons deposited in the detector. 

In this paper, we describe the search for neutrinos with energies above 10 6 GeV using data collected with the 
half-completed IceCube detector in 2008—2009. This analysis is sensitive to all three neutrino flavors. Compared to 
the previous EHE neutrino search described in Ref. (l3j , which used an earlier stage of the IceCube detector, the 
current analysis benefits from the enlarged instrumented volume and from improved agreement between simulated and 
observed event distributions. This article presents the improved strategies implemented since the previous analysis [13| . 



II. DATA SETS 

The analysis uses data collected from April 6, 2008 through May 20, 2009. At the time of data collection, the IceCube 
detector consisted of 2400 Digital Optical Modules (DOMs) on 40 vertical strings. The volume of the detector was 
roughly 0.5 km 3 with the detector center located at a depth of 1948 m below the ice surface. The DOMs consist 
of a 25 cm photomultiplier tube (PM T'l jl6| with data acquisition and calibration electronics, data compression, 
communications, and control hardware [17[ ■ The trigger setting was unchanged from the previous analysis (l3j ] . 

The analysis was optimized on simulated data with most of the experimental data kept blind. A 10% subset of 
the experimental data was used for examinations of the Monte Carlo simulations and detector response. This subset 
comprised 35.8 days of detector livetime distributed randomly throughout the data collection period, and was not used 
once the analysis was fully defined. The use of statistically independent final sample conservatively ensures avoidance 
of possible analysis bias due to tuning a Monte Carlo simulation using an experimental subset. The selection criteria 
were then applied to the complementary 90% of the experimental data, comprising 333.5 days of livetime. 

The primary background in this analysis is muon bundles made up of large numbers of muons produced by high 
energy cosmic-ray interactions in the atmosphere. This background was simulated with the CORSIKA air-shower sim- 
ulation package version 6.720 [l8| with the SIBYLL 2.1 [lj| and QGSJET-ii [2(| hadronic interaction models, without 
prompt muons from the heavy meson decays. Cosmic-ray interactions assuming pure proton and iron primary com- 
positions in the energy region between 10 6 and 10 10 GeV were simulated. Background contributions from primary 
cosmic-ray energies beyond 10 10 GeV were estimated by extrapolation of the simulated sample up to the GZK cut- 
off energy of ~5 x 10 10 GeV. EHE neutrino signal events in energies between 10 5 and 10 11 GeV from several flux 
models UHl, Hi], [22[ were simulated using the JULieT package [23[ . 



III. EVENT SELECTION 

The amount of energy deposited in the form of Cherenkov photons by the neutrino-induced charged particles in 
the detector is highly correlated with the energy of the particles [l3[ . An EHE neutrino interaction occurring inside 
or close to the IceCube detector would stand out against the background of cosmic-ray induced muons due to the 
much higher light deposition. The total number of photo-electrons (NPE) recorded in an event was used as the main 
distinctive feature to separate signal from background. 
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FIG. 1. Event observables in the quality bright sample that are used for the final selection criteria. Distributions of (a) NPE 
and (b) cosine of the reconstructed zenith angle for shallow events, and (c) NPE and (d) AtLN-E for deep events in a livetime 
of 333.5 days. The black circles represent experimental data and the solid and dashed lines are CORSIKA-SIBYLL with iron 
and proton primaries, respectively. The expected signal distributions from simulations of the GZK 1 model (sum of all three 
neutrino flavors) are shown as long-dashed histograms. Systematic uncertainties are not included. 



A. On-line sample 

The number of photo-electrons (p.e.) recorded by an individual DOM was derived by integrating the pedestal 
subtracted waveforms. Each DOM has two waveform digitizers, that simultaneously capture p.e. signals with differing 
dynamic ranges and time windows [17]. The event total NPE was then obtained by summing the number of p.e. 
detected by each DOM. PMT saturation effects and the sizes of the time windows limit the NPE estimation at high 
light levels. The initial NPE calculation was performed online at the South Pole. For this analysis we consider 
only events with NPE on i; nc > 630. The event rate of this "on-line bright sample" was ~1.4 Hz. At this level, the 
background rate exceeded the expected signal rate by > O(10 7 ). 



B. Off-line sample 

For the following data selection step, the NPE values were re-calculated after eliminating photon signals from low 
energy muons accidentally coincident in a 20 fis time window of a large NPE event. These low energy muons leave a 
faint light, typically with an NPE < 9. The light deposition of the coinciding low energy muon was, in most cases, 
spatially and temporally separated from the main bright p.e. cluster. While the few coincident photons have very 
small impact on the NPE calculation, they can disturb the geometrical reconstruction of the particle tracks later on 
in the analysis. Contributions from coincident low energy muons were eliminated by removing p.e. signals that were 
temporally separated from the time of the highest light deposition associated with the main high NPE event. The 
recording time of a p.e. signal in the i t h DOM, t Wyi , was defined as the time at which 10% of the total charge had 
been captured. The time of the highest light deposition was defined as the time tLN of the DOM which captured 
the largest p.e. signal in the event. This time, tLN , was typically associated with the time of closest approach of the 
charged particle tracks to any DOM in the detector. For the off-line NPE calculation and track reconstruction, those 
p.e. signals which occurred outside the time window [—4.4 /is, 6.4 fjs] around the tLN were excluded. The "off-line 
bright sample" selects events with NPE > 3.2 x 10 3 and the number of hit DOMs (NDOM) > 200; here and below 
NPE and NDOM are obtained after the t^jv time window cleaning. These NPE and NDOM thresholds reduced the 
background rate by two orders of magnitude while keeping ~70 % of the cosmogenic neutrino-induced events remain. 
The remaining backgrounds are bundles containing many hundreds of muons, with an estimated cosmic-ray energy 
above 10 7 GeV. 



C. Quality cut 

Apart from NPE, the particle direction and the depth distribution of the detected Chcrcnkov photons are distinctive 
event features that separate the EHE neutrino signal from the atmospheric muon background. Due to the energy 
dependence of the neutrino interaction cross section, most of the EHE neutrino signal is expected from directions 
close to the horizon. As a result of the depth dependence of the optical properties of the Polar ice, the largest photon 
signals are often detected in the deepest part of the detector, where the ice is most transparent 24]. On the other 
hand, the background atmospheric muons enter the detection volume from above and lose a substantial fraction of 
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their energy during propagation through the detector. Therefore, the time and depth coordinates, z, of the detected 
Cherenkov photons, measured relative to the detector center, show negative correlation for background. The largest 
photon signals from these background muons are expected at shallow depths near the top of the detector. Exceptions 
are inclined atmospheric muon bundles that pass outside the instrumented volume with the point of closest approach 
in the deep, clear ice at the bottom of the detector, or individual muons that deposit most of their energy in an 
isolated catastrophic energy loss in the deep ice after having passed through the top part of the detector. Track 
reconstructions often fail to identify such atmospheric muon events as downward-going tracks, when most of the light 
deposition occurs in the deep part of the detector. Therefore, a track reconstruction is applied only to those events in 
which the DOM with the largest signal is located at z > —300 m ("shallow events"). The negative z value indicates 
the vertical distance below the center of the IceCube detector. For events with the largest photon signal at z < 
—300 m ("deep events"), further event selection criteria rely on timing instead of directional information. 

For the shallow events, the particle directions are reconstructed with the LineFit algorithm [l3j]. Since the majority 
of the EHE neutrino induced events is close to the horizon [23| while the directions of the background muon bundles 
are mostly vertical, it is important to minimize the number of background tracks that are misreconstructed as 
horizontal. In order to reject the mis-reconstructed background events, another simple one-dimensional reconstruction 
is introduced. The distribution of average depth of p.c. as a function of timing arc fitted by a linear function, 
z(tio) — Co + Szt • ^io- The fit parameter, S z t, is a measure for the speed at which the light signal propagates 
in z-direction, and hence for the inclination of the tracks. For vertically downward-going relativistic particles, the 
quantity S z t/c takes values ~ — 1, where c is the vacuum speed of light, whereas close to horizontal tracks yield values 
Szt/c ~ 0. The shallow "quality bright sample" requires an additional condition of (S zt /c + cos6>) > —0.4 where 9 
is the reconstructed zenith angle from the LineFit. This condition excludes events for which the one-dimensional fit 
suggests a significantly more vertical downward-going geometry than the LineFit. Both signal and background are 
reduced by less than ^2% by this criterion. Figure [T] shows the distributions of NPE (panel (a)) and cos 9 (panel (b)) 
for experimental data, background and signal simulations in the quality bright sample. The distributions of CORSIKA- 
SIBYLL with an iron primary composition show a reasonable agreement with experimental data while the total event 
rates are 50% overestimated by simulation. The zenith angle reconstruction resolution of the shallow quality bright 
sample is ^1.4° RMS for muon bundle background and >^2.5° for signal. This is because the signal experiences 
more stochastic energy losses along with hadronic cascades at its interaction vertices. 

The deep bright events (Zln < —300 m) are mostly events that traverse the bottom edge of IceCube or are uncon- 
tained events that propagate or cascade below the detector. The inclination of these events tends to be reconstructed 
more horizontally than the true direction. The agreement between the simulation and experimental distributions 
improves with increasing NPE threshold values for these events. Events with NPE < 10 4 are discarded from the 
deep quality bright sample in order to achieve a reasonable agreement between experimental data and simulations. 
Since the majority of the EHE neutrino-induced events have NPE > 10 4 , the effect on the signal efficiency by this 
requirement is minimal. A fraction of 96% of background is rejected by the cut, while 91% of signal is retained. The 
panel (c) in Fig. Q] shows the NPE distributions from the deep quality bright sample. 



D. Final selection 

The final event selection is chosen in order to minimize the model discovery factor (MDF = /iids /^signal) [25] in 
the region of the phase space where a better signal to background ratio (S/B) is expected, where [i\<\ & is the least 
number of events to claim signal discovery at 5<r significance and N s i gna \ is the number of neutrinos expected from the 
GZK 1 [3] model flux. For the shallow events, high S/B is obtained in the region near the horizontal reconstructed 
direction as shown in Fig. [1] (b). For the deep events, instead of reconstructing the inclination of events, we utilize the 
time interval, AtLN-E, between the earliest detected photon in an event and tLN to obtain the best S/B subsample. 
The vertical atmospheric muon bundle events with the largest p.e. near the bottom of IceCube are often associated 
with a small number of p.e. in the shallow detector region much earlier than tLN- This contrasts to the EHE neutrino 
signal events. The main contributions to a detectable EHE signal in IceCube come from neutrino-induced horizontal 
muons and taus [23|. These produce the largest p.e. si gna ls shortly after the first recorded photo-electrons. Contained 
cascade-like events induced by neutrino interactions |26| inside the IceCube detector volume also exhibit a similar 
trend. Figure [11(d) shows the distributions in the deep quality bright sample. The best S/B is achieved in the bin 
AtLN-E ~ ns. The high rate in the experimental data for AtLN-E > 3600 ns is due to random noise in the DOMs 
and remaining coincident muons that were underestimated by the simulations. The slightly higher rate for the data 
in the bin AtLN-E ~ ns may reflect the fact that the ice is cleaner than what was simulated in the deep region. 
Figure presents the event distributions in the planes of cos 9 vs. NPE for the shallow events and AtLN-E vs. NPE 
for the deep events. Optimization is performed by differentiating the NPE threshold numbers in the region cos 9 < 
0.3 or AtLN-E < 0.5 /is for shallow and deep quality bright sample respectively. The NPE threshold of the other 
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TABLE I. Number of events passing cuts at various selection levels with 333.5 days detector livetime. The signal rates 
correspond to simulations of the GZK 1 model [3j. The errors of the on-line, off-line and quality bright samples are statistical 
only. Systematic uncertainties in the expected event rates at the final selection level are given as asymmetric error intervals 
after the statistical error. 



Samples 


Experimental 


Background MC (sibyll, iron) 


Signal MC (GZK 1) 


On-line 


3.7 x 10' 


(3.8 ±0.1) x 10'' 


1.8 ±0.007 


Off-line 


3.3 x 10 5 


(4.8 ±0.2) x 10 5 


1.2 ±0.006 




Shallow Deep 


Shallow Deep 


Shallow Deep 


Quality 


2.9 x 10 5 1.9 x W 


(4.4 ± 0.2) x 10 b (1.7 ± 0.2) x 10 J 


0.76±0.005 0.43±0.004 


Final 





0.076±0.012±°;°^ Q.032±0-QlQ±g;^ 


0.39±0.004t°;°g 0.18±0.002±°;°^ 



region (cos 9 > 0.3 or AtLN-E > 0.5 fis) is conservatively determined such that the number of background events 
above the threshold is less than 10~ 4 of the full livetime for each bin of cos 9 with width 0.2 or 1 /is for AtLN-E- 
This improves the detection sensitivity without sacrificing discovery potential. The solid lines in Fig. [2] are the final 
level selection criteria determined from the background (CORSIKA-SIBYLL, iron) and signal (GZK 1 [3]) Monte Carlo 
simulations following a blind analysis strategy. The minimum NPE threshold value is 2.5 x 10 4 . Events with NPE 
above the threshold value in each bin are considered to be signal event candidates. No events above the threshold 
are found in the 10% subset of the experimental sample. Monte Carlo simulations indicate that a cosmic-ray primary 
energy of at least ~2 xlO 9 GeV is required for a muon bundle to be selected as the final sample. Table Q] summarizes 
the number of events retained in each level of analysis. 




FIG. 2. Event number distributions of the shallow (upper panels) and deep (lower panels) quality bright sample in 333.5 days 
are shown for the background (left panels) and signal (right panels) simulations. The signal distributions are from GZK 1 
model [3] adding all three flavors of neutrinos. The background distributions are from CORSIKA-SIBYLL with iron primaries. 
The series of thick lines in each panel indicate the final sample selection criteria. 



IV. THE SYSTEM ATICS 

Table Hi] summarizes the sources of statistical and systematics errors in signal and background. The systematic 
uncertainties are assumed to have a flat distribution and are summed in quadrature separately for background and 
signal. 

The dominant source of systematic uncertainty in the signal event rate is the relationship between the measured 
NPE and the energy of the charged particles. The uncertainty is estimated by calibrating the absolute sensitivity of 
the DOMs in the laboratory and by calibrating the in-situ sensitivity using light sources co-deployed with the DOMs 
in the ice. The estimation by the latter method involves systematic errors in the simulation of the photon propagation 
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TABLE II. List of the statistical and systematic errors for signal (top) and background (bottom) simulations. The uncertainties 
for the signal are listed relative to the rate estimated for GZK 1 3]. The uncertainties in the signal rates vary with assumed 
signal spectra. The uncertainties in the background rate are estimated with CORSIKA-SIBYLL assuming iron composition. 



Sources 



Signal rate (%) 



Statistical error 

NPE 

Noise 

Neutrino cross section 
Photo-nuclear interaction 
LPM effect 



±0.8 
+3.9 / -7.2 
-1.8 
±9.0 
±10.0 
±1.0 



Total: ±0.8(stat.) ±14.0 -11.7(sys. 



Sources 



Background rate (%) 



Statistical error 

NPE 

Noise 

Cosmic-ray composition 
Hadronic interaction model 
Coincident events 



±17.0 
±37.1 / -46.7 
-2.2 
-83.9 
±36.1 
±31.2 



Total: ±17.0(stat.) ±60.4 -96.0(sys. 



in the ice. The uncertainty associated with possible underestimation in the DOM's random noise is estimated by 
adding artificial random pho to-electrons into 10% of the simulated events. The other uncertainties attributed to 
the neutrino interactions [27j and their daughters' interactions in the ice are similarly estimated as in the previous 
analysis [l3j . 

The dominant source of systematic uncertainty in the background event rates arises from the uncertainty in the 
primary cosmic-ray composition at the relevant energies (>10 7 GeV) and the hadronic interaction model used in the 
simulation of the air showers. The systematic uncertainty associated with the cosmic-ray composition is evaluated 
by considering two extreme cases of atmospheric muon simulations with either pure iron or pure proton primary 
compositions. Similarly, the uncertainty due to the hadronic interaction model is evaluated using atmospheric muon 
simulations with two different high energy hadronic interaction models: SIBYLL 2.1 and QGSJET-ii. Systematic 
uncertainties associated with the NPE measurement and the possible DOM noise rate underestimation are determined 
in the same manner as for signal events. The background contribution from possible prompt muons created in decays 
of charmed mesons is negligible. There is also uncertainty due to statistical limitations of the simulated coincident 
muon sample at NPE > 10 4 . This error is estimated by extrapolating distributions of statistically richer lower NPE 
coincident simulation events to the final selection region. Possible coincident events in the final sample are also 
estimated by the temporally and geometrically separated p.e. signals from the main p.e. cluster in each event. This 
coincident event check suggested that one of the two upward-going reconstructed events in Fig.[T]-(b) at cos 9 — —0.38 
was due to coincident muons. The other upward-going event (cos 9 = —0.83) was possibly neutrino-induced, while 
the NPE values of both events were approximately 4300 p.e., a factor of six less than the final threshold value. 



V. RESULTS 



No events in the blinded 90% experimental data pass all the selection criteria. This is consistent with the expected 
background level of 0.11 ± O^Iq'^q events in a livetime of 333.5 days. The passing rates for experimental and 
simulated events at each selection level are listed in Table |TJ 

The quasi-differential model-independent 90% CL limit on neutrino fluxes [33[ normalized by energy decade is 
shown in Fig. [3] assuming full mixing in the standard neutrino flavor oscillation scenario. In the limit calculation, 
the energy decade averaged effective area is used and the contribution from the Glashow resonance [29| is neglected. 
Incorporating the statistical and systematic uncertainties, the background is expected to be found with a uniform 
prior probability between and 0.19. This uncertainties are included in the final limit using a method outlined in [30j . 
This estimation together with the null result in the experimental sample gives the Feldman- Cousins 90% CL event 
upper limit [3l| of 2.35 events. For cosmic neutrinos with an E~ 2 energy spectrum, this implies an integral flux limit 
of E 2 4> < 3.6 x 10 -8 GeVcm~ 2 sec _1 sr _1 with the central 90% of the E~ 2 signal found in the energy interval 2.0 x 10 6 
— 6.3 x 10 9 GeV. This result is the first constraint of neutrino fluxes below the Waxman-Bahcall flux bound [32| in 
this energy region. 
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FIG. 3. All flavor neutrino flux differential limit and E 2 spectrum integrated limit from the 2008-2009 IceCube EHE analysis 
(red solid lines). The systematic errors are included. Various model predictions (assuming primary protons) are shown for 
comparison: GZK 1 ((m, Z max ) = (4,4)) @], GZK 2 Q, GZK 3 (fi A = 0.0), GZK 4 (fi A = 0.7) @, GZK 5 (maximal) and 
GZK 6 (the best fit, incorporating the Fermi-LAT bound) [6j. The gray dashed horizontal line indicates the Waxman-Bahcall 
flux bound with cosmological evolution ,321]. Model fluxes are summed over all neutrino flavors, assuming standard neutrino 
oscillations. The model independent differential upper limits by other experiments are also shown for Auger (PAO) 0, RICE [|, 
ANITA [ll[, and the previous IceCube result (IC22) [l3j]. Limits from other experiments are converted to the all flavor limit 
assuming standard neutrino oscillation and a 90% quasi-differential limit when necessary. The integral flux limit on a pure E~ 2 
spectrum is shown for AMANDA-II [121 ]. For reference, the estimated integrated limit for three years of observation with the 
full IceCube detector with the same analysis strategy is denoted as a gray solid line. 



TABLE III. Expected numbers of events in 333.5 days from several cosmogenic neutrino models and top-down models. The 
confidence interval for exclusion by this observations is also listed where appropriate. The cosmogenic neutrino models (GZK 1- 
6) assume the cosmic- ray primaries to be protons and different spectral indices/cutoff energies at sources as well as different 
cosmological evolution parameters and extension in redshift for the sources. Representative models with moderate (GZK 3, 4, 
6), moderately strong (GZK 1) and strong (GZK 2, 5) source evolution parameters are listed here. 



Models 


Event rate 


C.L. % 


GZK 1 [3] 


0.57 




GZK 2 [4J 


0.91 


53.4 


GZK 3 (fi A = 0.0) [5] 


0.29 




GZK 4 (fi A = 0.7) [5] 


0.47 




GZK 5 (maximal) [6] 


0.89 


52.8 


GZK 6 (the best fit) [6] 


0.43 




Top-down 1 (SUSY) [22J 


1.0 


55.7 


Top-down 2 (no-SUSY) [22] 


5.7 


99.6 


Z-burst [21] 


1.2 


66.4 


WB bound (with evolution) [32] 


4.5 




WB bound (without evolution) [32J 


1.0 





VI. DISCUSSIONS AND SUMMARY 



We analyzed the 2008-09 data sample collected by the 40-string IceCube detector to search for extremely-high energy 
neutrinos with energies exceeding 10 6 GeV. The differential and integral limits obtained are significantly improved 
relative to our previous result [13]. This is due to both the increased instrumented volume and improvements of the 
Monte Carlo simulations. The improved agreement between experimental and simulated data allowed a loosening of 
the NPE threshold in the data selection, thereby lowering the energy threshold of the analysis and improving the 
selection efficiency for high energy signal events that occurred outside the instrumented volume. This can also be 
seen in the corresponding neutrino effective area at the final selection shown in Fig. 2J Compared to the previous 
search [l3j], the effective area is a factor of 6 and 3.3 increased at 3 x 10 7 GeV and 10 9 GeV, respectively. The full solid 
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FIG. 4. Solid angle averaged neutrino effective area for four declination bands as well as that of the full solid angle average for 
(a) v e + v e , (b) Vy. + Pfi, and (c) v T + v T , assuming equal flux of neutrinos and anti-neutrinos. The lower right plot shows the 
final level signal event distributions for 333.5 days with the GZK 6 model spectra [f| for each neutrino flavor. 



angle averaged 3 flavor (assuming v e :v ^:v T — 1:1:1) neutrino effective area reaches 300 m 2 at 10 8 GeV and 2100 m 2 at 
10 9 GeV. The 90% CL differential limit at 10 9 GeV is a factor of ^4 higher than the fluxes predicted by the models 
GZK 2 and 5, and a factor of ^8 higher than the flux predicted by the models GZK 1, 4, 6, all of which assume 
primary protons. This suggests that the IceCube EHE neutrino search will reach these flux levels in the near future 
since the event rate is roughly proportional to the fiducial volume (see Fig. [3]), and the current analysis used only the 
half-instrumented IceCube detector configuration. Further improvements in sensitivity would enable IceCube to act 
as a a probe of the primary cosmic-ray composition at GZK energies [28j . 



Figure |4] indicates that a large part of the EHE neutrino signal are expected from the zenith angle region between 
60° and 90°. Upward-going EHE neutrinos are absorbed in the Earth. The propagation length of secondary muons 
and taus is greater than the distance between the surface and the IceCube fiducial volume. Thus, the inclined particles 
that reach the IceCube detector are created in the Earth. For v e , the event signatures are produced nearly at the 
neutrino interaction points and the current analysis is sensitive to all downward-going geometries. The peaked features 
in Fig. 4 (a) and (d) at E v<t ^6.3 PeV is due to the Glashow resonance [29]. Expected signal energy distributions of 
GZK 6 at the final selection level are shown in the lower right panel in Fig. SJ The peak energy of the expected signal 
after all selection criteria is at ~7.0 x 10 s GeV. Significant contributions from all neutrino flavors are observed. In the 
GZK 6 model, 13% of the signal are from v e , 45% are from and 42% are from v T . Through-going tracks (muons 
and taus) constitute 60% of the signal rate and the rest are neutrino interactions that create cascade-like events near 
and inside the detector volume. Table |Hl] gives the event rates for several model fluxes of cosmogenic neutrinos, 
top-down scenarios, and a pure E~ 2 power- law neutrino spectrum normalized to the Waxman-Bahcall flux bounds for 
reference. We expect 0.3 to 0.9 cosmogenic neutrino events in 333.5 days, assuming moderate to strong cosmological 
source evolution models. The half-instrumented IceCube detector is already capable of constraining those models 
with relatively high neutrino fluxes. The IceCube sensitivity to cosmological EHE neutrinos continues to grow. 
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Erratum: Constraints on the Extremely-high Energy Cosmic Neutrino Flux with the 
IceCube 2008-2009 Data 

An error has been found in the presentation of the neutrino effective areas in Figure 4 (a)-(c) of the original paper, 
which led to an overestimation of the values by a factor of six. Corrected neutrino effective areas are shown in Fig. [5] 
All other results reported in the paper, including the upper limit, were not affected and hence remain unchanged. 




FIG. 5. Solid angle averaged neutrino effective area for four declination bands as well as that of the full solid angle for (a) v e + v e , 
(b) + and (c) v T + v r , assuming equal flux of neutrinos and anti- neutrinos. 



